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Abstract 

We calculate one-loop corrections to the decays of the next-to-lightest neutralino %!] 
into the lightest neutralino Xi an d two leptons; this includes diagrams where a real 
photon is emitted. In cases where two-body decays X2 ~^ l^P —> Xi^~^ + are kine- 
matically allowed, we calculate these decays both with and without the single-pole 
approximation, and find consistent results. For example, for the minimal supergrav- 
ity parameter set SPSla, the integrated partial widths (the branching ratios) for 
X2 —> X?^ + (I = are enhanced by about 15.5 (13.6) percent by the one-loop 
corrections. We also study a scenario where \\ cannot undergo two-body decays, 
and find corrections to these branching ratios of about 13.6 percent. Moreover, we 
study the dilepton invariant mass (M;+^-) distribution, whose endpoint is often used 
in analyses that aim to reconstruct (differences of) supersymmetric particle masses 
at the LHC. The shape of this distribution is altered significantly by the emission 
of hard photons. For example, for the SPSla parameter set the peak of the Mal- 
distribution is shifted by several GeV when these contributions are included. 
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1 Introduction 



Supersymmetry (SUSY) [I] is one of the best motivated extensions of the Standard Model 
(SM) of particle physics. If SUSY exists at the electroweak scale, experiments at future 
high energy colliders should be able to discover the superpartners of known particles, and to 
study their properties [21 El- From the precise measurement of the masses, production cross 
sections and decay branching ratios of these superpartners, the fundamental parameters 
of the underlying SUSY models can be determined. This will help us to reconstruct the 
SUSY breaking mechanism. 

In the Minimal Supersymmetric Standard Model (MSSM) pQ with conserved i?-parity, 
the lightest supersymmetric particle (LSP), which in many scenarios is the lightest neu- 
tralino Xi, appears at the end of the decay chain of each supersymmetric particle. The LSP 
escapes the detector, giving the characteristic SUSY signature of missing energy. While 
this helps to suppress backgrounds from SM processes, it also makes the measurement of 
supersymmetric particle masses more difficult. This is true in particular at hadron colliders 
like the LHC, where the total energy in a given partonic collision is not known. 

At the LHC, the total SUSY production cross section is expected to be dominated by 
the production of gluinos and squarks, which decay into lighter charginos or neutralinos. 
Of particular interest are decay chains leading to the next-to-lightest neutralino Xt X2 * n 
turn can always undergo the three-body decays X2 ~^ Xiffi at least for light SM fermions 
/. Depending on the neutralino, sfermion and Higgs boson masses, the two-body decays 
X2 ~^ ff ~ > Xiff and/or X2 ~^ Xi^/cj) — > Xiff ma y also be open, where stands for one 
of the three neutral Higgs bosons of the MSSM; of course, Higgs intermediate states will 
contribute negligibly if / = e or \x. These leptonic final states are of particular interest, 
since they can be identified relatively easily even at the LHC. Moreover, the dilepton in- 
variant mass distribution can be measured accurately. In particular, the endpoint of this 
distribution is used in several analyses that aim to reconstruct (differences of) supersym- 
metric particle masses [21 B]. Under favorable circumstances it has been shown that this 
endpoint can be measured to an accuracy of 0.1% at the LHC [2]. In order to match this 
accuracy, at least one-loop corrections to x\ decays have to be included. 

Turning to the planned e + e~ linear collider ILC, X1X2 production is often the first 
process that is kinematically accessible [5] (other than Xi pair production, which leads to an 
invisible final state). The detailed analysis of X2 decays can then yield information about 
heavier supersymmetric particles. Under favorable circumstances, (9(10 4 ) X2 ~* Xil + l~ 
decays may be observed at the ILC, again making the inclusion of quantum corrections 
mandatory to match the experimental precision. In this paper we present a complete 
calculation of these corrections in the MSSM. 

The general MSSM has more than one hundred unknown free parameters. Therefore, 
it is not practicable to scan over the entire parameter space. Instead, several "benchmark 
scenarios" have been suggested p\, which are meant to illustrate characteristic features 
of various scenarios of SUSY breaking. Among those, the so-called SPSla parameter set, 
which has been defined in the framework of the mSUGRA scenario [L, has been studied 
particularly widely. It gives rise to a particle spectrum where many states are accessible 
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both at the LHC and at a 500 GeV ILC |4j. The masses of the relevant neutralinos 
and sleptons at this benchmark point are listed in Table ^ Note in particular that the 
two-body decays X2 ~ * lfl T — > x?Z~Z + are kinematically allowed; here 1% stands for the 
lighter one of the two charged sleptons of a given flavor. No other two-body decay mode 
is open. Moreover, squarks are so heavy that non-leptonic X2 decays can be neglected in 
this scenario. Note that Xi is mostly a U(1)y gaugino (bino), while X2 is dominated by 
its neutral SU(2) gaugino (wino) component; this is typical for most scenarios where the 
gaugino mass unification relation holds pQ. Leptonic two-body decays of X2 have been 
investigated at tree- level in Ref. three-body decays of X2 have been also studied at 
tree- level in Refs. j7J |S] . 



particle 


xl 


x\ 


e% (Ai) 


e 2 (£2) 




T~2 




v T 


mass [GeV] 


176.6 


96.2 


142.7 


202.3 


133.0 


206.3 


186.0 


185.1 



Table 1: Masses of the relevant neutralinos and sleptons for parameter set SPSla [0]. 

In this paper, we calculate leptonic X2 decays at one-loop level. Cases where X2 has 
two-body decays x\ ~^ ^ T ~^ xV- ^ + are treated both completely and in a single-pole 
approximation. In the complete calculation one has to employ complex slepton masses 
in the relevant propagators and one-loop integrals. The single-pole approximation in this 
case is performed by treating X2 decays as the production and decay of the sleptons l\. We 
compare the results of both methods, and find good agreement for the SPSla parameter 
set. We also analyze a scenario where X2 on by has three-body decays. In addition to 
calculating the integrated partial widths, we study the differential decay width of X2 as 
a function of the dilepton invariant mass. If X2 can undergo two-body decay, the shape 
of this distribution is essentially only affected by the emission of real photons; as well 
known, these contributions have to be added to the one-loop corrections to cancel infrared 
divergences. If x° only undergoes three-body decays, the shape of this distribution is also 
altered by the virtual corrections. In order to obtain the total decay width of X2 an d 
hence the branching ratios of its leptonic decays, the invisible decays x° ~~ * Xi u i^i an d the 
hadronic decays X2 ~~ > XiQQ are a l so calculated. 

The paper is organized as follows. In Sec. El we briefly summarize the renormalization 
of those sectors of the MSSM which are relevant for the decays of x\- The calculation of 
the tree-level decay widths is outlined in Sec. El Sec. E| discusses how to calculate these 
decays at one-loop level, including the emission of real photons. The complete one-loop 
calculation and the one-loop calculation in the single-pole approximation are presented in 
Sees. 14.11 and 14. 2\ respectively. In Sec. El the total decay width of X2 an d the branching 
ratios of the leptonic decays are studied. Some numerical results are given in Sec. El We 
conclude our work in Sec. [3 
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2 Renormalization of the MSSM 



In order to calculate the higher-order corrections, one must renormalize the parameters and 
the fields of the MSSM. Several approaches for the renormalization of the MSSM have been 
developed [3 El El El El El El- Here we employ on-shell renormalization following 
the strategy of Refs. [T2"l E] , This renormalization scheme is convenient for our purposes, 
since it ensures that the relevant supersymmetric particle masses are (almost) the same 
at one-loop level as at tree level; in particular, the endpoint of the Mu distribution is the 
same in both cases. We assume here that all relevant parameters are real quantities; this 
amounts to the assumption that the soft-SUSY-breaking terms conserve CP. 

2.1 Renormalization of the Chargino/Neutralino sector 

Loop corrections to the masses and mixing angles of charginos and neutralinos were first 
discussed in Ref. The independent SUSY parameters in the chargino/neutralino 

mass matrices are the electroweak gaugino mass parameters Mi, M 2 , and the supersym- 
metric Higgs mass parameter fi. These mass matrices also depend on the masses of the 
electroweak W and Z bosons, on the weak mixing angle $w, and on the ratio of vacuum 
expectation values (VEVs) tan/3; all these parameters are renormalized independently 
from the chargino/neutralino sector, as outlined below. In order to obtain finite S- matrix 
elements and Green's functions for chargino fields, we introduce a counterterm for the 
chargino mass matrix X, as well as field renormalization constants for the physical (mass 
eigenstate) four component (Dirac) chargino fields xf(i = 1,2) [T2] : 

X — ► X + 5X, (1) 

"Rxt — (fa + \ (SZ%^ u R xj , (2) 

where ul,r — (l-F7s)/2. Each element of 5X is the counterterm for the corresponding entry 
in X; in particular, its diagonal entries are the counterterms 5M 2 , 5fi. As for the fermionic 
fields of the SM, we need to introduce independent field renormalization constants for the 
left- and right-handed components of xt ■ These constants 5Z L and 8Z R are general 2x2 
matrices. 

Similarly to the chargino case, we introduce renormalization constants for the neutralino 
mass matrix Y and for the physical four-component (Majorana) neutralino fields x°(z = 
1,2,3,4) [T2]: 

Y 

^LX°i 
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Y + 5Y 



ulXj 



(3) 



(4) 



Here the elements of 6Y are the counterterms for the corresponding entries in Y; in par- 
ticular, the diagonal 2x2 blocks contain the counterterms 8M\, 5M 2 , and 5fi. The field 
renormalization constant SZ° is a general complex 4x4 matrix. Note that the Majorana 
condition x® = (x°) C implies that the left- and right-handed components of Xi do not 
renormalize independently, as shown in (J3J). 

In the on-shell renormalization scheme for the charginos/neutralinos ^2] the coun- 
terterms 5M 2 , 5/j,, and SMi are determined by requiring that the masses of xt j xt anc ^ 
Xi, which are defined as the poles of the corresponding propagators, are the same as at 
tree-level. We have slightly modified this prescription, keeping m^o, m^o and m-± fixed, 
since m^o is obviously more important for our analysis than Tn-±. The diagonal entries 
of the field renormalization constants are fixed by the condition that the corresponding 
renormalized propagator has unit residue. Furthermore, the renormalized one-particle 
irreducible two-point functions should be diagonal for on-shell external particles, which 
fixes the off-diagonal entries of the field renormalization constants. We note here that in 
this scheme the masses of the heavier neutralinos X3 4 an d lighter chargino xt do differ 
from their input (tree-level) values after one-loop corrections have been included. How- 
ever, these shifts violate the electroweak SU(2) symmetry, and are therefore usually quite 
small at least for gaugino-like states. 1 In most of mSUGRA parameter space (as well as in 
many other scenarios) the gaugino-like states are lighter than the higgsino-like ones. More- 
over, gaugino-like states are usually produced more copiously in the decays of squarks and 
gluinos. 

2.2 Renormalization of the Sfermion sector 

In general the superpartners Jl, fa of the fermions J'l, fa mix to form the sfermion mass 
eigenstates f s (s = 1,2). The MSSM does not contain right-handed neutrino superfields, 
hence there is no L — R mixing in the sneutrino sector. We assume that sfermions of 
different flavors do not mix. We renormalize the sfermion mass matrices Mj and the 

sfermion fields f s (s = 1,2) via 



The elements of the matrices 8Mj are the counterterms for the corresponding entries in Mj. 
The field renormalization constants 5Zj are general 2x2 matrices. For the sneutrinos, their 
masses, their counterterms, and the field renormalization constants are simple numbers 
rather than matrices. 

We follow Ref. ^3] and renormalize the sfermion sector via the on-shell scheme. For 
every generation of the squarks, the independent parameters are the soft-breaking mass 
parameters M~ L = Mj = M| i? M| H , Mj and the scalar trilinear coupling parameters A u , 

1 In the presence of strong L — B mixing in the stop sector, the masses of the higgsino-like neutralinos 
can still be shifted by several GeV |12l ITT) . 




Mi + SMr , 



(5) 




(6) 
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Ad. In order to fix their counterterms, one can renormalize two up-squarks and the lighter 
down-type squark via the on-shell renormalization scheme. It requires that the renor- 
malized masses of the two up-type squarks and the lighter down-type squark are equal 
to their physical (input) masses, and that the renormalized two-point function is diago- 
nal for on-shell external particles. These on-shell conditions determine the counterterms 
5M~ L , 5M~~ R , SMj as well as 5A U , 5Ad and the off-diagonal entries of the field renormal- 
ization constant, under the assumption SZg 12 = SZg 21 . The independent parameters for 
the sleptons are M~ , M~ and Ai. In analogy to the squarks, Ref. ^H] renormalized the 
sleptons via imposing the on-shell renormlaization conditions on the sneutrino and the 
lighter charged slepton. In a slight deviation from Ref. ^3] we fix both charged slepton 
masses at their tree-level values. In general there will therefore be a shift of the mass of the 
sneutrino when one-loop corrections are included. However, since this shift again vanishes 
in the limit of exact electroweak gauge symmetry it is numerically very small [T3J . Sim- 
ilarly to the chargino/neutralino case, the diagonal entries of SZj and the sneutrino field 
renormalization constant SZ^ are fixed by the requirement that the corresponding propa- 
gator has unit residue. Besides the soft-breaking sfermion mass parameters and the scalar 
trilinear coupling parameters, the sfermion mass matrices Mj also depend on fi, whose 
renormalization was discussed above, as well as on tan/3, m>z, 0w, the electric charge e and 
the charged fermion masses mf, whose renormalization will be discussed below. 



2.3 Renormalization of the neutral Higgs sector 

The renormalization of the Higgs sector in the CP- violating MSSM has been described in 
Ref. |14j : here we limit ourself to the neutral Higgs sector of the CP-conserving MSSM, 
using a mixture of on-shell and DR renormalization. 

The independent parameters in the Higgs sector are chosen to be the tadpoles T h o, T H o 
of the physical CP-even scalars h° and H°, which vanish at tree-level, the mass of the 
physical neutral CP-odd Higgs boson m A , and the ratio of VEVs tan/? introduced above. 
In addition the counterterms from the renormalization of the weak gauge boson sector, 
described below, enter here. 

In the neutral CP-odd Higgs boson sector, the mass matrix M x o and the fields A , G° 
are renormalized via 

M x o -> M x o + 5M x o , (7) 

A°\ ( 1 + \5Z AA \5Z AG \(A° 

G° J ^ \ \5Z GA 1 + \5Z GG J\G° 

Similarly, the neutral CP-even Higgs boson sector is renormalized as follows: 



M o -> M^o + dM^o, (9) 

h° \ ( 1 + \5Z hh \5Z hH \fh° 
H° J ^ \ \5Z Hh 1 + \5Z HH J\H° 



(10) 
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Mfjp is the mass matrix of the CP-even Higgs bosons and the matrices 5M x o } dM^o contain 
the counterterms 5T h o,5T H o, 5m\, and Stan/3. 

The tadpole counterterms are fixed by the requirement that the renormalized tadpoles 
vanish. The counterterm 5m A is determined by on-shell renormalization of the neutral 
CP-odd Higgs boson A . In this paper, we fix the field renormalization constants in the 
Higgs sector as well as S tan j3 in the DR scheme, which means that the counterterms only 
contain UV-divergent parts (plus some process-independent numerical constants). 2 In case 
of tan j3, this implies 

= ~n JlmZ A o z (m 2 A )] dtv . (11) 

tan p 2mz sin p cos p 

Other ways to renormalize tan/3 are discussed in Refs. [TK | FP U l27Ij. 



2.4 Renormalization of the SM-like sector 

The final piece of the Lagrangian we have to renormalize contains terms also present in 
the SM. Here we follow Ref. j2l]. The relevant parameters are the electric charge e, the 
charged fermion masses m/, and the masses of the W, Z bosons. They are renormalized as 
follows: 

e -> (l + 5Z e )e, (12) 
nif — > rrif + Srrif , (13) 

m w,z m w,z + Sm w,z- ( 14 ) 

The wave function renormalization of the fermion and neutral vector boson fields is de- 
scribed by 

ft - • \^!j)fj- (15) 

f? - fa + (16) 

(17) 



2 

Z\ ( 1 + Uz 7:7l Uz ZA \ ( z 



A ) \ \5Z AZ 1 + \5Z AA J \ A 



The renormalization constants above are again fixed by the on-shell conditions [21]. The 
on-shell definition of the weak mixing angle O w (sw = sin^vi/, Cw — cos^vf) is 



2 

2 _ ! m W 

o 



~ ± 2 



2 In the numerical examples discussed below, Higgs exchange contributions are negligible. However, 
they will be significant if ^2 ~~ * ^°Xi decays are open, and/or at high tan/3, where the Yukawa couplings 
to charged leptons and charge 1/3 quarks are enhanced. 
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Hence its counterterm is directly related to the counterterms of the gauge boson masses 

Ss w 



s w 



1 r 2 

2 s 2 



5m 2 w 



m 



5m 2 z 
mi 



(19) 



This completes our discussion of the renormalization conditions. We are now ready to 
discuss the calculation of the x° decay width. 
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3 Tree-level calculations for X2 — > x{l I 

The Born Feynman diagrams for \% — * X?^ + (^ = e ii l i T ) are displayed in Fig. ^ The 
propagators of the diagrams (a) and (b) have the structure as 

(20) 



k 2 



m 



2 ' 



where k and m\ denote the 4-momentum of the propagator and the slepton mass, re- 
spectively. If the two-body decays x® Xi^ ^ + are kinematically allowed, i.e. 
the sleptons l\ can be on shell at some points in the phase space, a finite width of l\ is 
necessary. It arises from the imaginary part of the slepton self-energy. A finite width is 
introduced via Dyson summation, 



k 2 — m 2 k 2 — ni- 
ii h 



it(k 2 



k 2 — m 2 

h 



+ 



k 2 -m 2 + E(k 2 ) 



(21) 



where S(fc 2 ) is the renormalized l\ self-energy. 

A gauge invariant matrix element is obtained by a Laurent expansion around the com- 
plex pole [2H]; in on-shell renormalization 



1 



k 2 - rrij + S(F 



1 - 



ReS(fc 2 ) 



mi 



k 2 



(22) 



mi 




Figure 1: The Born Feynman diagrams for x% — > Xil~l + {1 = e ii J -i T )- s = 1,2 labels 
the slepton mass eigenstates, denotes the MSSM neutral Higgs boson h°, H°, A , and 
the neutral Goldstone boson G° which appears only together with the Z boson in using a 
non-unitary gauge. Since the Yukawa coupling (f)l~l + is proportional to the lepton mass, 
the Higgs intermediate states are neglected when I = e and \i. 
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were m 2 v denotes the position of the complex pole in (J2TJ). It is obtained as the solution of 

m 2 p - m? + t{m 2 p ) = . (23) 

For the tree-level amplitude the complex pole ml is calculated at one-loop level. Its 
explicit expression is 

m 2 p = m{ - tmjTf^ , (24) 

where we have employed on-shell renormalization as in Sec. 12.21 r~ ree is the tree-level decay 
width of l\ and m^T-" 56 is the imaginary part of the slepton self-energy £(m~ ). The first 
factor in (J22j) is nothing but the Breit-Wigner propagator. Since the second term in the 
parentheses in (J22*|) is at one-loop level, we do not need it in the tree-level calculation. 
Therefore, the gauge invariant tree-level amplitude M tree for the decays x° ~~ * Xi^ + can 
be written as 

M tiee = * 2 1 9 + B(k 2 ), (25) 

k 2 - m 2 + imj T\ ree v ; ' v ; 

where K^of± ZT and ^"±|o/± re P resen t the xfyi^ an< ^ ^i Xi^ vertices, respectively, and B(k 2 ) 
denotes the non-resonant part of the matrix element, i.e. the matrix element of diagrams 
(a) and (b) for s = 2 and diagrams (c) and (d) in Fig. [T] 

The non-resonant part is much smaller than the resonant one (diagrams (a) and (b) 
for s — 1 in Fig. ^), hence it can be neglected approximately. We can then compute the 
relevant partial widths in the single-pole approximation, where the decays x® ~^ x\l~l + 
are treated as the production and decay of the sleptons /i, 

IXX2- XirZ + )t«e ^ r(*° - lfF) tree Br(lf - Xl^tree , (26) 

where the branching ratio of the decay if — > Xi^ 1S defined by 

Br(lf - Xi^ree = ^ ^ • ^ 

h 

The feature that the single-pole approximation reproduces the X2 partial width can be seen 
from the identity 

dk 2 , = L^r • (28) 



ptree 



' k 2 — m~ + iT\ Tee mr 
h h h 

If r~ rec <C m^ the integral in (|28|) will be dominated by the regions of k 2 close to m~ , i.e. 
only a narrow range of k 2 will contribute significantly. Moreover, the squared V~±%. ± is 

proportional to the if — > Xi^ partial width; together with the factor l/r~ ree from (|28j) 

this reproduces the factor Br(lf — > x?^) m 

Here we concentrate on scenarios where only the lighter charged sleptons l\ can be 
produced in x® decays; scenarios where sneutrinos and/or heavier charged sleptons can 
also be produced in two-body decays of X2 can be treated analogously. 
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4 One-loop calculations for X2 ~ * Xp ^ + 

The single-pole approximation can also be used at the one-loop level; however, we will first 
describe the complete calculation. 



4.1 Complete one-loop calculation 
4.1.1 Virtual corrections 

In general the virtual one-loop corrections to three-body decays can be classified as self- 
energy contributions, vertex contributions and box contributions. The first two classes 
are UV finite only after adding the contributions from the counterterms that originate 
from the renormalization of the MSSM, as discussed in Sec. 121 the box diagrams are by 
themselves UV finite. The MSSM Feynman rules, as well as the resulting counterterms, 
are implemented in the FeynArts package of computer programs [21], which allows an 
automated generation of the Feynman diagrams. The matrix element and the one-loop 
integrals are calculated with the help of the packages FormCalc and LoopTools [25] . 
respectively. 

Similarly to the tree-level case, diagrams with a slepton Zi propagator have singularities 
when l\ can be on shell. We remove the singularities by introducing a finite width of l\ as 
in (j21J) . Following the strategy in Sec. El one can obtain a gauge invariant matrix element 
at one-loop level. In order to obtain 0(a) accuracy near the li resonance, one needs to 
calculate the complex pole m 2 to two- loop level 

m l = m { ~ im h T i~ l00P > ( 29 ) 
where we have applied the on-shell renormalization scheme at two-loop level, and ri~ loop 

tl 

denotes the one- loop- level width of l\. Then the gauge invariant matrix element at one- loop 
level can be written as 

A(k 2 ) 

M tree + M virt = V oop + C(k 2 ) , (30) 

k z — m~ + im-j 1 - 
h 11 h 

where C(k 2 ) denotes the non-resonant part of the matrix element, the residue A(k 2 ) can 
be expressed as 

^PM fc2 )^S^ 2 ) + ^(^S^ 2 ) . ( 31 ) 

where V}^±°^ and V^}^ represent the renormalized one-loop corrections to the x^lfl^ 
and liXi^ vertices, respectively. 




9 



Figure 2: Examples for the virtual photonic corrections 



Feynman diagrams like those shown in Fig. El also give singularities when the sleptons 
li are on shell. The left (vertex) diagram has an infrared (IR) divergence if a real slepton 
mass is used in kinematic configurations where the slepton can be on shell. The right (box) 
diagram has a divergence which can be understood as being due to re-scattering of the two 
charged leptons in the final state, which persists for large photon virtualities. One should 
therefore use complex slepton masses, 



k 2 -mf k 2 -m 2 + imr T l f loop " ' 

in the one-loop integrals from these diagrams. This gives rise to a large QED logarithm 
log ^m^/ri~ loop y Furthermore, the box diagram shown in Fig. El has the property that 

the virtual photon is attached to external on-shell charged particles. This results in IR 
divergences, which we regularized by introducing a fictitious photon mass A. The IR 
divergences cancel after we add contributions from real photon bremsstrahlung, which will 
be discussed in Sec. 14.1.21 The one-loop integrals with complex masses can be calculated 
with the help of LoopTools. The analytical expressions for scalar three-point and four- 
point functions with real arguments can be found in Refs. j^U I2H I2Z|- We generalized 
these to allow for complex arguments. Note that in our calculation the masses of the light 
leptons, i.e. m ; (I = e,/i), are neglected except when they appear in the one-loop integrals, 
while the r mass m T is kept everywhere. 



4.1.2 Real photon bremsstrahlung 

In order to cancel the IR divergences in the virtual corrections, we have to add contributions 
from real photon bremsstrahlung, which contain on-shell propagators of stable particles in 
the limit where the scalar product of the 4-momenta of the photon and the emitting charged 
particle vanishes. 

This always happens, regardless of the mass of the emitting particle, if the photon 
energy is very small. The "soft photon bremsstrahlung" contribution is defined via 
the condition E 7 < AE; here the cutoff parameter AE should be small compared to 
the relevant physical energy scale (e.g. the resolution of the experimental apparatus). 
This IR-divergent contribution is sufficient to cancel the IR divergences from the virtual 
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corrections. Since the energy of the emitted soft photon is by definition very small, this 
emission essentially does not change the momenta of the other final state particles. This 
contribution is therefore described as a convolution of the differential tree-level decay width 
with a universal factor. Explicit expressions can be found in Refs. |2T| . 

Real emission contributions with > AE are called "hard photon bremsstrahlung" . 
Altogether, 

■Tbrems Tsoft 

(A£)+r hard (A£). (33) 

The dependence on the largely arbitrary parameter AE cancels after summing soft and 
hard contributions, provided it is sufficiently small. In the limit of vanishing mass of the 
emitting particle, the hard photon bremsstrahlung contribution also contains a divergence, 
if the momenta of the photon and the emitting particle are collinear. Since there are 
no massless charged particles in Nature, this is not a real divergence; in our case, it is 
regularized by the masses of the leptons in the final state. However, since the lepton 
masses, i.e. m e and m M , are very small, it is very difficult to get stable numerical results 
from a direct numerical evaluation of hard photon bremsstrahlung, e.g. using Monte Carlo 
integration. 

This can be overcome by dividing hard photon bremsstrahlung into a collinear part, 
where the angle between the photon and the radiating particle is smaller than a very small 
angle A9, and the complementary non-collinear part, 

r hard (A£) = r coll (A£,A0)+r non _ coll (A£,A0). (34) 

The angular cutoff AO should be so small that the emission of photons emitted at angle 
< AO relative to the emitting lepton can be assumed not changing the direction of the 
3-momentum of this lepton. 

If we treat a charged lepton and a collinear photon inclusively, i.e. the momentum 
of collinear photon is added to that of the emitting lepton, analytically the differential 
contribution of the collinear emissions is written as the differential tree-level decay width 
multiplied by a universal function j2Hll2ni- This approach is for collinear-safe observables 
[2H] . If one adds the soft and collinear contributions to the virtual corrections, all singu- 
larities (lnm; and In A) cancel. This is in accordance with the Kinoshita-Lee-Nauenberg 
theorem jSU]. At the LHC the electron energy is determined calorimetrically: in this case 
a collinear photon would hit the same cell of the calorimeter as the lepton, so the two en- 
ergies cannot be disentangled. 3 Hence the electron observables are defined as collinear-safe 
observables in our calculation. 

We also consider non-collinear-safe observables [2H], where the lepton and its collinear 
photon are not treated inclusively. The contribution of the collinear photon bremsstrahlung 
cannot be calculated analytically. In this case the mass singularity lnm^ cannot be canceled 

3 There is a minor caveat to this statement. The experimental definition of an "electron" usually requires 
the existence of a charged track whose energy - more exactly, absolute three-momentum - should not be 
grossly different from the energy measured by the calorimeter. This requirement may remove a few events 
with very hard collinear photons. 
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in the differential width and hence becomes visible. Muon energies are generally measured 
through the curvature of their track in a magnetic field. This measures the energy (more 
exactly, the 3- momentum) of the muon after emitting the collinear photon (if any). In our 
calculation the muon observables are treated as non- collinear- safe observables. Finally, 
the contribution from the emission of hard non-collinear photons is calculated by using a 
multi-channel Monte Carlo approach |31j . 

The virtual photonic corrections by themselves are UV divergent 4 , hence one cannot 
meaningfully separate the QED corrections from the one-loop contributions by simply 
selecting diagrams which contain a photon. In the case of the light lepton (/ = e, fi) final 
states, following conventions of the Supersymmetric Parameter Analysis (SPA) [32], we 
can pick out and separate potentially large QED terms from the the sum of virtual and 
soft photon bremsstrahlung corrections, r virt + r so ft : 



virt 



Tsoft r -|- Premainder , (35) 



where V contains all the potentially large terms proportional to logm; or logAi?, while 
Premainder is IR and UV finite and free of such large QED logarithms. The "QED contri- 
butions" can then be defined as follows: 

TqED = T + Thard • (36) 

Note that Tqed defined in this way does not depend on the cutoff parameters AE and 
Ad. Moreover, terms proportional to logm; cancel between the two contributions in Tqed 
in (jSEJ) (specifically, between T and r co ii) in the integrated width and in the differential 
width for the collinear-safe observables. Using the definitions (j3*3j) and ()36|) . the complete 
one-loop contribution can be written as 

Pi— loop Ttree Tvirt Tbrems 

= r tree + r remam der + TqED • (37) 



One should perform the replacement (}32jl also in the real photon bremsstrahlung 
when the sleptons l\ can be on shell. In this case V contains the large QED logarithm 
log ^m^ri~ loop j , besides logm; and \ogAE. However, in the integrated partial width 
these terms nearly cancel after summing all contributions; more exactly, the total pre- 
factor of log {ra-^/T 1 ^ 100 ^ vanishes when ri~ loop _^ o, once one includes the fact that the 

squared ifxlP vertex is oc T^. 5 

When t~t + are the final states of \X decay, the r mass m T is kept everywhere. This 
mass is so large that a stable numerical result can be obtained from the hard pho- 
ton bremsstrahlung even for A6 — > 0, i.e. we do not need to divide the hard photon 

In principle one can define a renormalizable non-supersymmetric theory containing only leptons, slep- 
tons, neutralinos and photons. However, the counterterms computed in this theory would be different 
from those of the full MSSM. 

5 In the limit Tj i — > some kinematical distributions would become singular; for example, the distribu- 
tion in the invariant mass of the Xi ~ ^ systems would contain 6— functions. 
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bremsstrahlung contribution into collinear and non-collinear parts. Since we do not count 
terms oc log "large logarithm" , we follow a slightly different procedure to define the 

"QED part" of the correction. The virtual corrections contain photonic and non-photonic 
contributions, 

r V irt = r 7 irt + r™" 7 , (38) 

both of which are UV divergent, while the sum is finite (after including all counterterms) . 
The photonic virtual corrections can be split into an UV-finite part T and an UV-divergent 



part T 7 



UV-div' 

Tvirt = r + r 7 JV _ div . (39) 



Here r 7 JV _ div contains the terms that would be subtracted in an DR regularization of r 7 irt . 
After this rearrangement, the virtual corrections can be written as 

r — f _l_ F 7 _|_ pnon-7 

virt — 1 Tl UV-div 1 virt 

r "I - rrcmainder j (40) 

where r remainder = r 7 JV _ div + r™°™~ 7 as well as T are UV finite. The "QED corrections" are 
finally defined as 

Tqed = r + 

L brcms j 

(41) 

where Tbrcms stands for the contribution from all diagrams with real photon emission. By 
construction, Tqed is both UV and IR finite. 



4.2 One-loop calculation in the single-pole approximation 

If ~^ hi two-body decays are allowed and X2 does not have other two-body decay modes, 
at one-loop level, just like at tree level, the decays X2 ~ * Xil~l + can be approximately 
treated as production and decay of the sleptons l±, 

T(X 2 - X?^ + )i-ioop * T(x° 2 - itni-ioapBrtft - X?/ ± )i_ioop , (42) 

with 

Mh -> Xil )i-ioop = • (43) 

The virtual contributions of the production and decay of the sleptons l±, which now only 
contain vertex type corrections but no box diagrams, are again calculated with the help 
of the programs FeynArts, FormCalc and LoopTools. In order to obtain IR-finite 
results, the real photon bremsstrahlung is added, which is again separated into an IR- 
divergent soft part and an IR-finite hard part. For the light lepton final states I — e, /i, 
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the division of the hard photon bremsstrahlung contribution into a collinear part, which 
can be calculated analytically, and a non-collinear part, which is calculated numerically, 
proceeds along the lines described in Sec. 4.1.2. As discussed in Sec. 14. 1| the UV-divergent 
photonic contributions cannot be treated separately as "QED corrections" . We define the 
"QED corrections" in the same way as in the complete calculation. One finally arrives at 
a total one-loop contribution which is independent of the cutoff parameters. 



5 Total decay width of jjQ an d the branching ratios of 
the decays jjQ ~ > X?^ + 

As discussed in Sec. [TJ the next-to-lightest neutralino X2 can decay into the LSP x% an d 
two fermions //. The leptonic final states are important because they can be identified 
at the LHC Moreover, the endpoint of the dilepton invariant mass distribution is used to 
determine the mass relations of supersymmetric particles. The invisible x\ decay modes, 
i.e. X2 ~ * Xi u i^h do not effect the dilepton invariant mass distribution. But they contribute 
to the total width of yQ- Since it is very difficult to identify quarks at the LHC, the hadronic 
decays X2 ~ > XiQQ are l ess interesting than leptonic decays. In order to obtain the total 
decay width of x°> these hadronic decays must be calculated. The total decay width of X2 
can be written as 

r x§ = E M - ri+ *S) + r (*2 - iwxS)] + E r (x° - MX?) • ( 44 ) 

l=e,n,r q=u,d,c,s,b 

Here we assume that the decay X2 ~^ Xitt * s n °t kinematically allowed. The branching 
ratios of the leptonic decays X2 ~^ Xil + l~ are defined as 

Br(x° 2 - X\l + l-) = ^ ~Z *° ll+n ■ (45) 

The invisible decays x9> ~ > x\ v \Pi are calculated at tree and one-loop level. At tree 
level these decays can proceed through the exchange of sneutrinos or Z bosons; cf. Fig. 1. 
Here we focus on the case where the decays x\ ~^ x\ v i v l are pure three-body decays. 
These decays are calculated similarly to the calculations for the leptonic three-body decays 
X2 ~ * Xil~l + - Since none of the external particles carries electric charge in the decays 
X2 ~^ Xi u l u h there are no corrections involving real or virtual photons, and hence no IR 
divergences. Therefore, there are also no QED corrections in these decays. This makes the 
calculation of the partial width into neutrinos considerably simpler than for decays into 
charged leptons. 

The hadronic decays of X2 are calculated in order to obtain the total width of The 
Born Feynman diagrams for the decays X2 ~^ XiQQ (<? 7^ t) are similar to those of Fig. 1 
with intermediate squarks instead of sleptons. Here we only consider the case where the 
hadronic decays X2 ~^ XiQQ are pure three-body decays. Since the SUSY-QCD corrections 
are not considered in our calculations, these decays can be treated in the same way as 
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X® — > x^ ^ + - The virtual photonic corrections of the hadronic decays (the diagrams are 
similar to the ones of the leptonic decays, i.e. Fig. |2J) are IR divergent. The contributions 
of the real photon bremsstrahlung are necessary for the cancellation of the IR divergences. 
We neglect the light quark masses, i.e. m q (q = u, d, s), except when they appear in the one- 
loop integrals. In analogy to Sec. 14.1.21 the contributions of the real photon bremsstrahlung 
are also splitted into an IR-divergent soft part and an IR-finite hard part. For the light 
quark final states, we separate the hard photon bremsstrahlung into a collinear part and a 
non-collinear part in order to obtain stable numerical results. Since quarks are detected as 
jets, which contain many photons, quark energies are always collinear-safe. As presented 
in Sec. I4.1.2[ the soft and collinear contributions are calculated analytically. We treat the 
decays with heavy quark final states in the same way as t~t + final states. The QED 
corrections are defined in the same way as in Sec. 14.1.21 since the photonic contributions 
are UV divergent and cannot be treated separately. 

6 Numerical results and discussion 

We are now ready to present numerical results of our calculation. We present results both 
for a scenario where x[> can undergo two-body decays, and for a scenario where no two-body 
decays of x° are possible. Furthermore, we discuss both the integrated partial widths and 
branching ratios of X21 an d the distribution of the invariant mass; this distribution is 
of great interest for future experiments, as discussed in the Introduction. Since we always 
assume equal masses for selectrons and smuons and the light lepton mass mi (I = e, fi) is 
neglected except when it appears in the one-loop integrals, the integrated partial widths for 
the e + e~Xi and n + ji^Xi fi na l states are identical. As discussed in Sec. I4.1.2| the dilepton 
invariant mass M e + e - is defined as collinear-safe observable, i.e. we add the momentum 
of a collinear photon to that of the emitting electron, since it is difficult to separate their 
energies at the LHC The energies of a muon and its collinear photon can be disentangled 
easily at the LHC, hence the dilepton invariant mass M M + M - is defined as non-collinear-safe 
observable, i.e. the momentum of a collinear photon is not added to that of its emitter 
muon. In this case the large logarithm In m M can not cancel in the distribution, so the mass 
effect can be seen in the dilepton invariant mass distribution. One will obtain identical 
M e + e - and M M + M - distributions if both of them are defined as collinear-safe observables. In 
order to see the differences of the two treatments (adding and not adding the momentum 
of a collinear photon to the emitting lepton), we also show the comparison of dilepton 
invariant mass M^^- and M e + e - distributions. 

6.1 Numerical results for the SPSla parameter set 

We first present results for the SPSla benchmark scenario, as described in Table 1 in the 
Introduction. The dilepton invariant mass M e + e - distribution from the complete calcula- 
tion is shown in Fig. El In the left frame we show not only the tree-level and total one-loop 
predictions, but also the separate QED and "remainder" corrections, see ()37|) . We see that 
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Figure 3: The dilepton invariant mass M e + e - distribution for the SPSla parameter set 
(complete calculation). 



the non-QED contributions are positive and quite large everywhere, whereas the QED con- 
tribution is large and negative near the endpoint of the distribution, but small elsewhere. 
In full three-body kinematics this endpoint is simply given by l3_bod y = m x% ~ m x?" 

However, for the SPSla parameter set X2 decays are dominated by contributions with on- 
shell l\ in the intermediate state. The endpoint for this two-body configuration is given by 



^eT e X -| 2 _ body =m,0, 



mz± 



m 2 - 

Xi_ 

2 

m + 



76.8 GeV, (46) 



where the numerical value holds for the SPSla scenario. Note that this is only 3.6 GeV 
below the endpoint of the three-body decays. At tree level, the M e + e - distribution peaks 
at the region which is a little below the endpoint of the two-body contribution. The right 
panel in Fig. |31 which shows a blow-up of the endpoint region, shows that the peak of 
this distribution is then moved about 4 GeV below the endpoint (j4*6^) once higher-order 
corrections are included. This is almost entirely due to contributions where a hard photon 
is emitted, which takes away energy from the e + e~ system. This change of the shape 
of the invariant mass distribution near the endpoint is important, since in (simulated) 
experiments one needs a fitting function describing this distribution in order to determine 
the location of the endpoint [33]. in Fig. El we used A6 = 1° in the definition of collinear 
photons. In a real experiment, even photons emitted at somewhat larger angles might be 
counted as contributing to the energy of the emitting electron. In this case the change of 
the shape of the M e + e - distribution will be somewhat smaller. 

In Fig. 0] we compare the numerical results of the complete calculation and the single- 
pole approximation at tree (left) and one-loop level (right). At tree level the M e + e - distri- 
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M e+e -(GeV) M e+e -(GeV) 

Figure 4: The dilepton invariant mass M e + e - distribution for the SPSla parameter set 
(comparison of the complete calculation with the approximate calculation). 



bution computed in the single-pole approximation has an exactly triangular shape, with a 
sharp edge at the endpoint ([46)1 . This edge is smeared out a bit in the complete tree-level 
calculation, which includes the full set of diagrams shown in Fig. 1. As noted above, this 
edge is also softened considerably once hard photon emission is included. The single-pole 
approximation therefore works even better in the one-loop calculation. However, this excel- 
lent agreement even for the differential decay width is partially accidental. The agreement 
would become somewhat worse if the endpoints in two- and three-body kinematics were 
further apart; this would happen if the mass of l\ was close to the mass of either or to 
that of Xi, since then one of the two square roots in (f4T)J) would become small. 

The comparison of the dilepton invariant mass M At + At - and M e + e - distributions is shown 
in Fig. E| In the upper frames we show the dilepton invariant mass Mi+i- (I = e, fi) 
distribution both at tree and one-loop level. Since the selectrons and smuons have equal 
masses and the light lepton mass mi (/ = e, /i) is neglected except when it appears in the 
one-loop integrals, their distributions are identical at tree level and different at one-loop 
level due to the different treatment of the collinear-photon radiation. From these figures 
one obtains that at one-loop level the mass effect is larger near the endpoint than in other 
regions and the peak of the M^+^- distribution is shifted to lower invariant-mass values in 
comparison with the M e + e - distribution. We also show the relative one-loop corrections 
in the lower frames in Fig. E] The relative one-loop corrections from the final state 

is smaller than that of the e + e~ final state in the upper invariant-mass region, while it is 
larger in the lower invariant-mass region. The main reason is that we add the momenta 
of collinear photons to that of emitting electrons, but we do not do this for the collinear- 
photon radiation from muons. Hence the invariant mass M^+^- is reduced in comparison 
with M e + e -. This leads to the shifting of events from the upper invariant-mass region to 
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Figure 5: The comparison of the dilepton invariant mass M At + M - and M e + e - distribution 
for the SPSla parameter set. 



the lower invariant-mass region. 

The corresponding results for the t + t _ x? final state are shown in Figs. HI and [7| Fig. EI 
shows that both the QED and, in particular, the non-QED corrections are smaller in 
magnitude than for light leptons. In case of the QED contribution this is essentially a mass 
effect. Our angular cutoff A6 defining the collinear region in (J34)) is so small that even 
non-collinear radiation off electrons or muons is still more likely than any hard radiation 
off t leptons; recall that we do not split hard radiation into collinear and non-collinear 
contributions for t + t^Xi final states. 

The reduction of the non-QED corrections is even more dramatic. They amount to 
about +20% for electrons and muons, but only to about +6% for tau leptons. This 
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M T+T - (GcV) M T+T - (GeV) 

Figure 6: The dilepton invariant mass M T + T - distribution for the SPSla parameter set 
(complete calculation). 

difference stems from the fact that l\ is a pure £77(2) singlet for I = e, //, since we neglect 
terms oc mi in the mass matrices of these sleptons. In contrast, ?i — tr mixing is quite 
significant, leading to a sizable SU(2) doublet component of f\. Therefore X2 decays into 
(real or virtual) l\ can only proceed through its small U(l)y gaugino (bino) component 
for I = e, fi, while the large SU(2) gaugino (neutral wino) component also contributes 
for I = t. Moreover, the Xi^T^i coupling, which is involved in the virtual corrections, 
only exists for I = r. The reason is that li is a pure SU(2) singlet for I — e, fi while 
f\ has a sizeable SU(2) doublet component, as explained above. In the limit of exact 
SUSY, these contributions involve the U{1)y and SU{2) gauge couplings, respectively, 
which renormalize (and run) quite differently. If there are significant differences between 
masses of supersymmetric particles, the differences between the true gauge couplings and 
these gaugino-lepton-slepton couplings also becomes significant [Slj; note that in the SPSla 
scenario, squarks (which contribute to various two-point functions) are about three times 
heavier than \X- Finally, the x^i 7 " vertex also receives non-negligible contributions which, 
again in the limit of exact SUSY, are proportional to the r Yukawa coupling [lj. As a 
result of the reduced non-QED corrections, the total correction is now negative, especially 
for large values of M T + T — . 

The endpoint region of the M T + T - distribution from the complete calculation is shown 
in the right panel of Fig. |U1 We see that the peak of the distribution is shifted downwards 
by about 2 GeV once higher-order corrections are included. A shift of this magnitude 
may be significant, even though the r + r~ invariant mass is in general difficult to measure 
accurately, due to the presence of v T (anti-)neutrinos in the r decay products, which escape 
detection. 

In Fig. [7| predictions from the complete calculation are compared to those from the 
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Figure 7: The dilepton invariant mass M T + T ~ distribution for the SPSla parameter set 
(comparison of the complete calculation with the approximate calculation). 

single-pole approximation. In this case we find almost perfect agreement even in the 
endpoint region, both at tree level and after including one-loop corrections. The reason is 
that for SPSla, m fl happens to be very close to ^om^o. Performing the replacement 
— > in ()46|) shows that the endpoints of the r + r~ distributions in two- and three- 
body kinematics practically coincide. 

The partial widths of the different X2 decay modes and the branching ratios of its visible 
leptonic decays are listed in Table where the numbers in the parentheses are obtained 
from the approximate calculations. We find — > Xi p i^i) ^ J2i^i.X2 ~* Xi^ + ^)- 
This is not surprising, since the charged lepton final state is accessible via on-shell l\ 



decay mode 


tree-level width(MeV), Br 


one loop-level width(MeV), Br 


e~e + Xi 


1.123 (1.122), 5.9% 


1.297 (1.294), 6.7% 


V~V + Xi 


1.123 (1.122), 5.9% 


1.297 (1.294), 6.7% 




16.870 (16.933), 88.0% 


16.595 (16.646), 86.2% 


VeVeXl 


0.012 


0.012 




0.012 


0.012 




0.013 


0.013 


qqxl (q^t) 


0.015 


0.015 


total width 


19.168 


19.241 



Table 2: Partial widths of different X2 decay modes and the branching ratios of its visible 
decays for the SPSla parameter set. The numbers in parentheses give the corresponding 
partial widths calculated in the single-pole approximation. 
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intermediate state, whereas for the neutrino final state all exchanged particles are off shell. 
Since squark masses are near 500 GeV in SPSla scenario, hadronic final states contribute 
even less than neutrinos do. 

From the results in Table El one concludes: 

• The main decay mode of x\ * s X2 ~^ T T+ Xv Its branching ratio is about 88.0% at 
tree-level, 86.2% at one-loop level. This mode dominates partly because of the lower 
mass of fi as compared to e\ (133.0 GeV vs 142.7 GeV). Even more important is that 
the X2^i r coupling is much stronger than the x\e\e coupling, which in turn is due to 
significant L — R mixing, which only exists in the f sector, as explained above. 

• The total decay width is enhanced by 0.4% when one- loop corrections are included. 
Such modest corrections are typical in the absence of large enhancement factors (e.g., 
large logarithms). This overall perturbative stability confirms that our choice of 
renormalization scheme, and of the electroweak input parameters listed in Appendix 
A, is indeed rather well suited for the task at hand. 6 

• One-loop corrections enhance the partial width and the branching ratio of X2 ~^ 
l~l + Xx = e '/ i ) decays by 15.5% and 13.6%, respectively. This results from the 
large size of the positive non-QED corrections depicted in Fig. 3. Much of these 
corrections can probably be absorbed into an appropriately defined running 
coupling. This is illustrated in Fig. |HJ which compares the Mj+j- distribution com- 
puted including only the "universal corrections" defined in Ref. |TT] (see also |34j ) 
with the tree-level and full one-loop results. We see that the residual non-universal 
corrections are relevant only close to the edge of the lepton pair distribution, where 
real photon emission is most important. Since this result is for a specific scenario, a 
more comprehensive analysis might be appropriate. 

• The single-pole approximation reproduces the integrated partial widths to about 
0.3% accuracy. This agreement is even better than in the M;+ z - distribution shown 
in Figs. 4 and 6. In fact, from (}2*%)) and the discussion at the end of Sec. 4.1.2 
one might expect better agreement for the integrated partial width than for (some) 
kinematical distributions. 

6.2 Numerical results for pure three-body decays 

We also investigated the effect of higher-order corrections on leptonic X2 decays for a 
scenario where x\ does not have any two-body decay modes. To that end we again use the 
SPSla parameter set, except that the soft SUSY-breaking parameters in the slepton mass 
matrix are set to 

m lh = 230 GeV , m~ lR = 183 GeV, I = e, /i, r . (47) 

6 Of course, the total width after the inclusion of one- loop corrections is scheme independent, up to 
unknown two-loop correction terms. However, the relative size of the one-loop corrections does depend on 
the chosen scheme. 
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Figure 8: Dilepton invariant mass M^^- from an approximate calculation with only uni- 
versal 1-loop contributions via effective couplings, in comparison with the tree-level and 
the complete 1-loop results. SPSla parameter set. 

The masses of the relevant neutralinos and sleptons in this modified SPSla parameter set 
are listed in Table El where one finds that X2 has to undergo a pure three-body decay. 
Therefore we do not have to introduce complex slepton masses in the one-loop functions. 
Apart from this simplification, the calculation is very similar to the "complete" calculation 
described in Sec. 4.1. 



particle 


xi 


X°i 


ei (fix) 


e 2 (A2) 


Tl 


T 2 


ui(l = e,/x,r) 


mass (GeV) 


176.6 


96.2 


187.9 


234.9 


182.3 


239.2 


221.0 



Table 3: Masses of the relevant neutralinos and sleptons for the modified SPSla 
parameter set. 

The dilepton invariant mass M e + e - and M T + T - distributions are shown in Figs. 03 and 
ITUl respectively. At tree level the M e + e - distribution shows a small peak near its upper 
endpoint from the exchange of nearly on-shell Z bosons. Since the QED and non-QED 
corrections are very small and negative in this region, this peak is less pronounced once 
one-loop corrections are included. This is of some significance, since the shape of this 
distribution can now be used to infer the strengths of various contributing diagrams, which 
in turn provides information on slepton masses and neutralino mixing ESI- Since f 



1 1 1 1 1 1 1 r 
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exchange is much enhanced relative to e exchange, one cannot see any contributions of 
Z-exchange even at tree level from the M T + T ~ distribution. Moreover we can observe 
that the invariant mass M e + e - and M T + T - distributions have a rather sharp edge at their 
endpoints. These edges are again softened by real photon emission, but remain quite 
distinct. This should facilitate the experimental determination of the endpoint, and hence 
the measurement of m^o — m^o . 

We compare the dilepton invariant mass M M + M - and M e + e - distributions in Fig. fill 
where the tree- and one-loop-level results, the blow-up of the endpoint region and the 
relative one-loop corrections are shown. From these figures one obtains that the shapes of 
the M M + M - and M e + e - distributions are identical at tree level and different at one-loop level 
due to the different treatment of collinear-photon radiations. In contrast to the numerical 
results from the SPSla parameter set (see Fig. EJ), the mass effect is small in Fig. [TTJ but 
it is still distinct, especially in the relative one-loop corrections. In the calculations for the 
invariant mass distribution, the momentum of a collinear photon is added to that of the 
emitting electron, but it is not added to that of the emitting muon. Hence the invariant 
mass M M + M - is reduced in comparison with M e + e -. It leads to the shifting of events from 
the upper invariant-mass region to the lower invariant-mass region. This effect can be seen 
in the lower frames in Fig. ^2 i- e - m the lower invariant-mass region the relative one-loop 
corrections of the final state is larger than that of e + e~ final state, while the inverse 

relation holds in the upper invariant-mass region. 

The decay width of different X2 decay modes and the branching ratios of its visible 
decays are shown in Table EJ The total X2 decay width is about 600 times smaller than for 
scenario SPSla. This is not surprising, since in this scenario X2 can only have three-body 
decays, while the two-body decays X2 ~~ * n ^ ~~ * Xi^ ^ + are kinematically allowed in the 
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Figure 9: The X2 decay width differential in the dilepton invariant mass M e + e - in the case 
of genuine three-body decay. 
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dM T+T - 




Figure 10: The x 2 decay width differential in the di— r invariant mass M T + T - in the case 
of genuine three-body decay 

scenario SPSla. In the modified SPSla scenario one-loop corrections increase the total x% 
decay width by a modest 1.2%. 

Turning to the various partial widths of leptonic decays, we notice that the f + f~Xi 
final state is still the largest decay mode of x\ (25.1% at tree level, 25.4% at one-loop level) 
since is smaller than the selectron and smuon masses and the large L — R mixing exists 
only in the f sector. Note that exchange of the SU(2) doublet sleptons now dominates for 
I = e, (j, since the size of the y^e^e coupling exceeds that of the \%^R e coupling by nearly a 
factor of 10 and all the sleptons are off shell in this scenaro. This dominance of cl exchange 
also explains why the e + e~Xi an d ^e^eXi final states now have quite similar partial widths: 
in the limit where x\ an d Xi are pure 577(2) and U(l)y gauginos, respectively, the product 
of couplings involved in and v e exchange is exactly the same (up to an overall sign). 

The pattern of one-loop corrections to leptonic decays is different from the original 
SPSla scenario. The partial width into electrons and muons is still enhanced by about 
14.3%. But now the invisible partial widths are also increased, diminishing the correction 
of the branching rations. In the original SPSla scenario, the one-loop partial widths 
of the invisible decays are almost unchanged in comparison with the tree-level ones, see 
Table El Note that we assumed three exactly degenerate sneutrinos here, unlike in the 
original SPSla scenario, where v T is slightly lighter than v e . In the modified scenario a 
tiny difference between the partial widths for v T v r x\ and v e v e x\ final states nevertheless 
results from one-loop corrections involving the r mass or Yukawa coupling (e.g. from the 
v and v two-point functions). 

The hadronic final states have very large partial decay widths and branching ratios: 
r hadronic = 15 -210 keV (52.9%), rj^°°£ ic = 14 - 765 keV (50.7%), though the squark masses 
are much larger than the slepton masses. Part of the reason is that the Z-exchange dia- 
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Figure 11: The comparison of the dilepton invariant mass M At + /J - 
in the case of a genuine three-body decay. 
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grams give larger contributions to hadronic final states than to leptonic ones. Moreover, 
the interference between Z and sfermion exchanges is large and positive for the hadronic 
final states, while it is also large but negative for the leptonic final states. This is the main 
reason why the hadronic decays of x\ obtain so large branching ratios. 



25 



decay mode 


tree-level width(keV), Br 


Hoop- level width(keV), Br 


e _ e + y ( i l 


1.270, 4.4% 


1.451, 5.0% 


A, 1 


1.270, 4.4% 


1.451, 5.0% 


t~t + y? 

/V. 1 


7.209, 25.1% 


7.383, 25.4% 


v P v P y^ 


1.273 


1.355 




1.273 


1.355 




1.273 


1.354 




2.480 


2.386 


ddxi 


3.330 


3.298 


CCXi 


2.475 


2.378 


ssxi 


3.330 


3.298 


bbxi 


3.595 


3.405 


total width 


28.778 


29.114 



Table 4: The decay width of different X2 decay modes and the branching ratios of its visible 
leptonic decays in the modified SPSla scenario. 

7 Summary and conclusions 

We have performed a complete one-loop calculation of the decays x® ~^ ^ ^ + Xi (I — 
e, fj,, r). The necessary renormalization is briefly described in Sec. 2. In most cases we 
used on-shell renormalization, which leaves the masses of the relevant neutralinos and 
sleptons (almost) unchanged. This is convenient for our purpose, since one important goal 
in the experimental study of leptonic X2 decays is the determination of (differences of) 
super symmetric particles masses from the dilepton invariant mass (Mj+j-) distribution. 

For the cases where the intermediate charged sleptons can be on shell, these decays were 
calculated both completely and in a single-pole approximation at one-loop level. In the 
complete calculation one has to employ complex slepton masses in the relevant propagators 
and one-loop integrals. The single-pole approximation in this case is performed in the way 
that the x° decays are treated as a sequence of two two-body decays. We checked that 
for the well-studied SPSla parameter set, this approximation reproduces the integrated 
partial widths to better than 0.5% accuracy even after one-loop corrections are included. 
For this parameter set we find a rather small one-loop correction to the total x\ decay 
width, but the branching ratio for the most easily detectable electron and muon final 
states are increased by about 13.6% at one-loop level. 

We also studied the effect of higher-order corrections on the Mj+j- distribution. If 
only one exchanged particle can be on-shell, as in the SPSla scenario, the shape of this 
distribution is altered only by real photon emission contributions, i.e. its peak is shifted by 
several GeV below the endpoint. This is very important since the shape of the distribution 
near the endpoint should be known if the endpoint is to be determined accurately from real 
data. In our calculation we define collinear photons as being emitted at an angle A9 < 1° 
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relative to the emitting lepton. Since the selectrons and smuons have equal masses and 
the light lepton mass mi (I = e, /i) is neglected except when it appears in the one-loop 
integrals, one will obtain identical distributions for M e + e - and M^+^- if the momentum of 
a collinear photon is added to that of the emitting lepton. The actual effect of the collinear- 
photon radiation depends on details of the measurement apparatus, and therefore has to 
be calculated anew for each experiment. We have focused on the LHC experiment in our 
calculation. At the LHC the electron energy is determined calorimetrically. In this case 
a collinear photon would hit the same cell of the calorimeter as the electron, so the two 
energies cannot be disentangled. Hence we add the momentum of a collinear photon to the 
one of the emitting electron in our calculation. Since muons pass through the calorimeter, 
where the photons are detected, and measured further outside in the muon detector at the 
LHC (their 3-momenta are measured through the curvature radius of their track in the 
magnetic field), the momentum of a collinear photon is not added to the one of its emitter 
muon in our calculation. In this case the mass effect can be seen in the dilepton invariant 
mass distribution. We find that the peak of the M e + e - distribution is moved downwards 
by about 4 GeV once the one-loop corrections are added. In comparison to the M e + e - 
distribution, the peak of the M M + M - distribution is shifted slightly to lower invariant-mass 
values at one-loop level. This is due to the different treatment of the collinear-photon 
radiation. 

We have also analyzed a modified SPSla scenario, with increased slepton masses, so 
that X2 can only undergo genuine three-body decays. In this case we again find a moderate, 
if slightly larger, correction to the total x\ width when one-loop corrections are considered, 
but the branching ratios for the electron and muon final states are still enhanced by about 
13.6% at the one-loop level. We have seen in Fig. [H] that for the simpler case / = e,/i 
the bulk of the non-QED correction to the partial width can be absorbed into new xV-il 
couplings, which are sensitive to the spectrum of sfermions. In the case of r final states, 
significant ft, — Tr mixing as well as the r Yukawa coupling have to be included in the 
analysis. We have not attempted to define such effective couplings and, perhaps, mixing 
angles here. 

In this modified SPSla scenario, the dilepton invariant mass distributions have a rather 
complicated shape, showing the contributions from Z exchange near the upper endpoints 
of the distributions. In this case the shape of these distributions is affected not only by 
real photon emission, which again leads to significant negative corrections for large M1+1- , 
but also by virtual corrections, which can e.g. differ for Z and slepton exchange diagrams. 
In this case the shape of the distribution away from the endpoint also carries information 
about slepton masses and neutralino mixing angles. Fitting tree-level distributions to real 
data might therefore give wrong results for these physical parameters. In this context 
a careful analysis of collinear radiation is also important, since differences in the energy 
measurements of electrons and muons could lead to spurious differences of fitted selectron 
and smuon masses. Here the collinear-photon radiations for electrons and muons are 
treated as discussed beforehand. One finds that the one-loop shapes of the M e + e - and 
M A1 + M - distributions are different, though the selectrons and smuons have equal masses in 
our calculations. 
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We conclude that higher-order corrections to leptonic x° decays can exceed the 10% 
level both in integrated partial widths and branching ratios, and in the shape of the dilepton 
invariant mass distribution. Attempts to absorb much of the large virtual corrections into 
effective running couplings might be rewarding. An accurate understanding of x° decays 
is of considerable importance, since this is one of the lightest visible particles that can be 
produced directly at future e + e~ colliders, and plays a prominent role in the analysis of 
cascade decays of gluinos and squarks at the LHC. 
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Appendix 



A Parameters 

For the numerical evaluation, the following values of the SM parameters are used: 

m e = 0.510999MeV, m M = 105.6584MeV , m T = 1.777GeV , 

m u = 53.8MeV, m c = 1.5GeV, m t = 175GeV, 

m d = 53.8MeV, m s = 150MeV, mf, = 4.7GeV, 

m w = 80.45GeV, m z = 91.1875GeV , 

a(0) = 1/137.0359895, G M = 1.1663910 x lCT 5 GeV~ 2 . 

The on-shell renormalization scheme requires a = a(0) for one-loop calculations. For the 
tree level expressions we instead use the effective coupling for the overall normalization, 



2 



a Gr = ^Mk!K (4g) 



7T 



We saw in Sec. 5 that this choice leads to good perturbative stability of the total x° decay 
width. 
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